Entropy of isolated quantum systems after a quench 
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A diagonal entropy, which depends only on the diagonal elements of the systems density matrix in the en- 
ergy representation, has been recently introduced as the proper definition of thermodynamic entropy in out-of- 
equilibrium quantum systems. We study this quantity after an interaction quench in lattice hard-core bosons 
and spinless fermions, and after a local chemical potential quench in a system of hard-core bosons in a superlat- 
tice potential. The former systems have a chaotic regime, where the diagonal entropy becomes equivalent to the 
equilibrium microcanonical entropy, coinciding with the onset of thermalization. The latter system is integrable. 
We show that its diagonal entropy is additive and different from the entropy of a generalized Gibbs ensemble, 
which has been introduced to account for the effects of conserved quantities at integrability. 
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The notion of entropy was first used by Clausius in the mid- 
19th century and was soon put in the context of statistical me- 
chanics by Boltzmann and Gibbs. Generalized to quantum 
mechanics by von Neumann in the 1930s and incorporated by 
probability theory by Shannon in the 1940s, entropy has man- 
ifested itself in different forms over the years. Despite the di- 
versity, the consensus is that any physical definition of entropy 
must conform with the postulates of thermodynamics QjJ, |2EI - 

An appropriate definition of entropy, suitable also for iso- 
lated quantum systems out of equilibrium, is fundamental for 
advances in nonequilibrium statistical mechanics and for a 
better understanding of recent experiments with quasi-isolated 
quantum many-body systems, such as those realized with ul- 
tracold atoms fl. von Neumann's entropy, defined as Sn — 
— Tr(plriyS), where p is the many-body density matrix (the 
Boltzmann constant here and throughout this Letter is set to 
unity), complies with the laws of thermodynamics when de- 
scribing isolated quantum systems in equilibrium and quan- 
tum systems interacting with an environment, but it becomes 
problematic when dealing with closed systems out of equilib- 
rium. Since in an isolated system Sn is conserved for any 
process, this entropy is not consistent with the second law of 
thermodynamics. This motivated the recent introduction of 
the diagonal (d) entropy |4|. which is given by 



Sd = - Pnn In(Pnn), 



(1) 



where p nn are the diagonal elements of the density matrix in 
the instantaneous energy basis. In equilibrium Sd coincides 
with the von Neumann's entropy. In addition, Sd was argued 
to satisfy the required properties of a thermodynamic entropy: 
it increases when a system in equilibrium is taken out of equi- 
librium, it is conserved for adiabatic processes, it is uniquely 
related to the energy distribution (and as such satisfies the fun- 
damental thermodynamic relation), and it is additive. 

More specifically, it was indicated in 14J that Sd should be 
equivalent to the equilibrium microcanonical entropy when 
the energy fluctuations are subextensive and the energy dis- 



tribution is not sparse, assumptions expected to hold in non- 
integrable systems. For integrable systems, the existence of 
a complete set of conserved quantities JHt] invalidates those 
assumptions and precludes thermalization in the usual sense. 
However, it has been shown that few-body observables after 
relaxation can still be described by a generalized Gibbs en- 
semble (GGE) which is a grand-canonical ensemble ac- 
counting for the conserved quantities ||7|] . 

Here, we study the d entropy in isolated quantum systems 
after a quench in both integrable and nonintegrable regimes . 
We consider two kinds of quenches in one dimension (ID): an 
interaction quench for hard-core bosons (HCBs) and spinless 
fermions, which have a nonintegrable regime [6||, and a local 
chemical potential quench for HCBs (or spinless fermions) 
with a superlattice potential, which are integrable. In the 
first case, as the system transitions to chaos, we show that 
the distribution function of energy becomes Gaussian-like and 
Sd approaches the thermodynamic entropy. This indicates 
that thermodynamically the system becomes indistinguishable 
from a thermal state. In the second case, Sd is shown to be ad- 
ditive and different from the entropy of the GGE. This differ- 
ence scales linearly with the system size, suggesting the exis- 
tence of additional correlations not captured by the GGE . 

Quench and entropies. We consider a particular initial 
state \ipim) which is an eigenstate of a certain initial Hamil- 
tonian. At time r = 0, the Hamiltonian is instantaneously 
changed (quenched) to a new one with eigenstates |^„) and 
eigenvalues E n . The initial state then evolves as \i>{r)) = 



J2n Cn 



,-iE„ 



where C n = (^ n \ipmi) and |C„| 2 cor- 



respond to the diagonal elements, p nn , of the density matrix, 
p(r) = |V(t))<V(t)|. 

For generic systems, with nondegenerate and incommensu- 
rate spectra, the expectation values of few-body observables 

(O) relax to the infinite time average (0(t)) = ^2 n p n nO nn , 
which depends o nly on the diagonal elements p nn and O nn = 
(*„|6|* n ) CSEill- Thus, the d entropy (HJ is the entropy of 
the diagonal ensemble. It resembles the Shannon entropy, but 
with no arbitrariness in the basis. For sudden quenches, Sd is 
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equivalent to Sn for the time averaged density matrix. The 
difference between Sd and thermodynamic entropies can help 
to quantify additional information contained in the diagonal 
part of the density matrix and not in the equilibrium ensem- 
ble. 

One may also write Sd as the sum of a smooth S s and a 
fluctuating S f part Sd = S s + Sf 0], where 

S s =J2PnnHv(E n )SE], (2) 

n 

Sf = ~y] prm ln[p nn rj(E n )5E]. (3) 

n 

Here r){E n ) is the density of states at energy E n : rj(E) = 
^2 n S(E — E n ) and SE 2 is the energy variance: SE 2 = 
J2 n Pnn(E-E ini ) 2 , where E ini = (V>ini|# IV'ini) is the expec- 
tation value of the quenched Hamiltonian with respect to the 
initial state. In the continuum limit, S s = j dEW(E)S m (E) 
and S f = - J dEW(E) ln[W(E)6E], where W(E) = 
J2 n PnnS(E — E n ) is the energy distribution. In S s , the 
microcanonical entropy, S m (E) = \n[?](E)SE], is the log- 
arithm of the total number of accessible states in the range 
of energy [E - 6E/2,E + SE/2]. If the system is large 
and finite-size effects become negligible, then up to nonex- 
tensive corrections, S m becomes equal to the canonical en- 
tropy, S c = -E n i Z ~ le ~ En/7 'ln(Z- 1 e- iJ "/ T )], where T 
is the temperature related to the energy of the system and 
Z = Y^ n e~ En / T is the partition function (see Ref. Il2ll ). 

When W(E) is narrow, so that SE is subextensive, S s be- 
comes equivalent to the equilibrium microcanonical entropy. 
Moreover, if in addition W(E) is a smooth function of en- 
ergy, then Sf is also subextensive. These features are ex- 
pected to be generic for the nonintegrable (chaotic) regime, 
where the eigenstates (away from the edges of the spectrum of 
systems with few-body interactions) become pseudo-random 
vectors SB. 

In the integrable limit, on the other hand, conserved quanti- 
ties reduce the number of eigenstates of the Hamiltonian that 
have a nonzero overlap with the initial state ITll , SO p nn 
becomes sparse and Sf extensive. In this case, both terms S s 
and Sf are expected to contribute to the d entropy. It then 
becomes appropriate to compare Sd with the entropy of the 
GGE introduced in Ref. J(J, which accounts for the integrals 
of motion. The many -body density matrix of the GGE is given 
by pgge = Z G ^ E e-^ x ^\ where Z CCE = Tr[ e -£ A ™H 
{J TO } is a complete set of conserved quantities, and A m are 
the Lagrange multipliers fixed by the initial conditions A m = 
ln[(l - (V'ini|i : m|V-'ini))/('0ini|i : m|V'ini>]- Since the GGE is a 
grand-canonical ensemble, which can suffer from large finite 
size effects for small systems, in addition to the entropy in the 
GGE, Sgge, we also compute the entropy in its canonical ver- 
sion (GCE) as the trace Sgce = Tr[/5 G cE ln(/3 G cE)]can where 
only eigenstates of the Hamiltonian with the same number of 
particles contribute to the trace. 

Chaotic systems. We consider periodic ID chains with 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 



hopping and interaction, with the following Hamiltonian 

L 



H 



n i+ 2 2 



(4) 



for hard-core bosons and similarly for spinless fermions (with 



( , and hj — > rtp, where standard notation 



T, and -> hi 
has been used (81]. L is the lattice size and we take the num 
ber of particles to be N = L/'i. We use full exact diagonal- 
ization to compute all eigenstates of the Hamiltonian, taking 
advantage of conservation of total momentum k due to trans- 
lational invariance. The initial states considered are eigen- 
states of Eq. with parameters ii n j, Vj n i, t', V belonging to 
the k = subspace. The final Hamiltonian (after the quench) 
has t = V = 1 and the same initial values of t' = V. The ini- 
tial states are selected such that their energies E m \ in the final 
quenched Hamiltonian are the closest to E at a chosen effec- 
tive temperature T, computed as E = Z^ 1 J2 n E n e~ E "/ T . 
When t' = V' = the system is integrable, while the addi- 
tion of NNN terms eventually induces the onset of chaos HI. 

Full exact diagonalization of the models above limits the 
system sizes that can be studied to a maximum of 8 particles 
in 24 lattice sites and thus prevents proper scaling studies of 
the entropies with increasing system size. This is left to the in- 
tegrable quenches where larger lattices can be explored. Here 
we compare Sd, S s , Sf, S m , and S c for the two largest system 
sizes available and for various Hamiltonian parameters as one 
departs from the integrable point. 

The main panels in Fig. Q] depict Sd and S s for systems 
with L = 24 at different effective temperatures as t', V in- 
creases. An agreement between Sd and S s can be seen as one 
approaches the chaotic limit, improving with temperature and 
system size [cf. insets in Fig. da) and[TJc)]. (By compar- 
ing the left and right panels, particle statistics does not seem 
to play much of a role.) Lower temperatures, for which Sd 
and S s are seen to depart, imply initial states whose energies 
are closer to the edge of the energy spectrum. For finite sys- 
tems, thermalization has been argued not to occur in those 
cases HI], and, from our results here, we expect that the idea 
of a thermodynamic description will break down if the tem- 
perature is sufficiently low. Increasing the system size is ex- 
pected to increase the region of temperatures over which a 
thermodynamic description will be valid. Figure[T]also shows 
that different initial states give slightly different quantitative 
results (top vs bottom panels), although the overall qualitative 
behavior is the same. 

The insets in Fig. [TJb) and Old), depict a comparison be- 
tween Sd and the equilibrium entropies in thermodynamic en- 
sembles whose energy has been chosen to be the same of the 
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FIG. 1: (Color online) Entropies vs t' — V . Left: bosons; right: 
fermions; top: quench from t\m = 0.5, Vm = 2.0; bottom: quench 
from tj n j = 2.0, Vi„i = 0.5. Filled symbols: d entropy |QJ; empty 
symbols: S s (3;OT = 1.5; □ T = 2.0; AT = 3.0. All panels: 
1/3-filling and L — 24; insets of panels (a) and (c) show Sd/S s for 
L = 24, thick (red) line, and L = 21, thin (black) line for T = 3.0. 
Solid lines in the insets of panels (b) and (d), from bottom to top: 
microcanonical entropy; canonical entropy S c for eigenstates with 
k — and the same parity as the initial state; S c for eigenstates with 
k = and both parities; and S c for all eigenstates with N — 8. 



initial state after the quench. Explicit results for the micro- 
canonical entropy with SE determined by the energy uncer- 
tainty are in surprisingly good agreement with those of Sd- 
Up to a nonextensive constant, the canonical entropy S c can 
also be written in the same form as S m © if we use the canon- 
ical width SE 2 = —dpE. Results for S c are shown for three 
different sets of eigenstates: (i) all the states in the N sec- 
tor, (ii) only the states in the N sector with k = 0, (iii) only 
the states in the TV sector with k = and the same parity as 
the initial state. The latter, as expected, is the closest to S m 
(also computed from eigenstates in the same symmetry sec- 
tor as iV'ini)) and Sd- In the thermodynamic limit, all three 
sets of eigenstates should produce the same leading contribu- 
tion to S c , but for finite systems it is necessary to take into 
account discrete symmetries in order to get an accurate ther- 
modynamic description of the equilibrium ensemble. 

The fact that Sd/S m — > 1 in the chaotic limit and that the 
agreement improves with system size provide an important in- 
dication that Sf is small and subextensive. Information con- 
tained in the fluctuations of the density matrix becomes neg- 
ligible in chaotic systems and only the smooth (measurable) 
part of the energy distribution contributes to the entropy of 
the system. Also, the close agreement between Sd and S m in 
the insets of Fig.QJb) andQJd) suggests that Sd is indeed the 
proper entropy to characterize isolated quantum systems after 
relaxation. Results for the energy distribution W(E) in Fig. [2] 
further support these findings. 

Figure [2] shows W(E) for HCBs for quenches in the in- 
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FIG. 2: (Color online) Normalized distribution function of energy. 
Bosonic system, L — 24, T = 3.0 and the values of t' = V' are 
indicated. Top panels: quench from tm = 0.5, Vi„i = 2.0; bot- 
tom panels: quench from tm = 2.0, Vm = 0.5. Solid smooth line: 
best Gaussian fit (v / 27ra) _1 e~( B ~ 6) ~' /(2a ~ ) for parameters a and b; 

-(E-E mi ) 2 /(2SE 2 ) 



dashed line: (V27r<5.E) e 

tegrable (left) and chaotic (right) domains. The sparsity 
of the density matrix in the integrable limit is reflected 
by large and well separated peaks, while for the noninte- 
grable case W(E) approaches a Gaussian shape similar to 
(V^kSE)-^-^- 11 '^ 2 ^ 26 ^, as shown with the fits. The 
shape of W(E) is determined by the product of the average 
weight of the components of the initial state and the den- 
sity of states. The latter is Gaussian and the first depends on 
the strength of the interactions that lead to chaos, it becomes 
Gaussian for large interactions [15]. A plot of p nn vs energy, 
on the other hand, does not capture so clearly the integrable- 
chaos transition lfl2ll . 

Integrable systems. We consider a ID HCB model with NN 
hopping and an external potential described by, 



L-l 



H S = -t H.c.) + A ]T cos ( y 



2*3 \ d 



b) bj . 



(5) 



This model is exactly solvable as it maps to spinless noninter- 
acting fermions (see e.g., Ref. lfl6ll ). The period P is taken 
to be P = 5, t = 1, and the amplitude A takes the values 
4, 8, 12, and 16. We study systems with L = 20, 25 .. .55 
at 1/5 filling. For the quench, we start with the ground state 
of © with A = and evolve the system with a superiattice 
(A 0) and vice-versa. Open boundary conditions are used 
in this case. 

We first study how the deviation of Sd from S s , as quanti- 
fied by Sf/Sd, scales with increasing lattice size for different 
quenches. As shown in Figs. [3] (a) and (b), Sf/Sd, does not 
decrease as L increases, rather, we find indications that Sf/Sd 
saturates to a finite value in the thermodynamic limit. Hence, 
for these systems Sd is not expected to be equivalent to the 
microcanonical entropy. 

In the lower panels of Fig. [3] we study the scaling of Sd with 
increasing system size for the same quenches. A clear linear 
behavior is seen, demonstrating that Sd is indeed additive. In 
these panels, we also show the microcanonical (with SE de- 
termined as for the interaction quenches fl7| ) and canonical 
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FIG. 3: (Color online) Entropy vs system size. Panel (a): from 
top to bottom, quench to — 4, 8, 12, 16; panel (b): quench 
from Ami — 4, 8, 12, 16, curves closely superpose. Lower pan- 
els: the quench type is indicated as (initial yl)-(final A). Symbols: 
Sd\ dashed lines: GCE-entropy (the closest to the d entropy in all 
cases studied); dotted lines: GE-entropy; dashed double-dotted line: 
canonical entropy; and dash-dotted line: microcanonical entropy. 

ensembles. The latter two can be seen to increase linearly with 
L and with a similar slope. These two entropies are clearly 
greater than Sd indicating that the diagonal ensemble in this 
case is highly constrained. Finally, we show results for the 
GGE and GCE entropies. They also increase linearly with 
system size and with a similar slope, showing that in the ther- 
modynamic limit their difference should be subextensive. In- 
terestingly, the slopes of the GGE and GCE are greater than 
the slope of the diagonal entropy. This suggests the existence 
of additional correlations not fully captured by the generalized 
ensemble. The diagonal entropy in this case is a clear observ- 
able independent measure of such correlations. This finding 
opens an important question as to which ensemble should be 
appropriate to characterize the thermodynamic properties of 
isolated integrable quantum systems after relaxation follow- 
ing a quench and for which observables these additional cor- 
relations are relevant. 

Summary. We presented a study of the diagonal entropy 
following quenches in integrable and nonintegrable isolated 
quantum systems. In the nonintegrable regime, we showed 
that Sd has the properties expected from an equilibrium mi- 
crocanonical entropy. In particular, the fact that Sd coincides 
with S m up to subextensive corrections and is thus determined 
only by the energy of the system implies that basic thermo- 
dynamic relations can be applied to nonintegrable isolated 
systems (see also discussion in Ref. J3l). In the integrable 
limit, we demonstrated that Sd is additive and smaller than 
the entropy of generalized ensembles (recently shown to prop- 
erly describe observables after relaxation following a quench). 
Our results open further questions as to how to characterize 
the thermodynamic properties of isolated integrable systems, 
and also motivate further studies for nonintegrable systems, in 



order to verify the scaling of Sd with system size and compare 
it to the one of the entropy in conventional statistical ensem- 
bles. 
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EQUILIBRIUM (SMOOTH) ENTROPY IN THE GAUSSIAN 
APPROXIMATION 



FUNDAMENTAL THERMODYNAMIC RELATION IN 
FINITE SIZE THERMALLY ISOLATED SYSTEMS 



The equilibrium thermodynamic entropy or the smooth part 
of the diagonal entropy for general non-equilibrium distribu- 
tions can be simplified if the energy fluctuations are small and 
the energy distribution is approximately Gaussian. In equilib- 
rium all ensembles (canonical, grand-canonical, microcanoni- 
cal, etc) assume that the density matrix is diagonal with prob- 
ability distribution being a smooth function of energy (with 
possible exceptions of the boundaries of the distribution like 
in the microcanonical ensemble). Therefore, for the equilib- 
rium ensembles, in the expression of the entropy one can sub- 
stitute the summation over the discrete energy levels by an in- 
tegration over the continuous energy spectrum. Starting from 
the von Neumann's entropy (which coincides with the diago- 
nal entropy in equilibrium), as in Eqs. (2) and (3), we obtain: 



dEW(E)S m {E) - J dEW{E)\n{W{E)5E), 

(6) 

where S m (E) = \n(r](E)SE) and rj(E) is the many-body 
density of states. Using the Gaussian shape of the energy dis- 
tribution: 

(E - Ef 



W(E) 



1 

/2n8E 



exp 



2SE 2 



(7) 



and expanding the entropy up to the terms vanishing in the 
thermodynamic limit, we find 



S m (E) +ln(V27r) + 



1 - 



SE 2 \ 



(8) 



where 5E 2 = — dpE is the equilibrium canonical width of 
the energy distribution. From this expression it is easy to 
check that the equilibrium entropy is indeed maximized for 
the canonical width SE = 5E C . Besides, this expression 
suggests that in the microcanonical definition of the entropy 
S m , the most natural width of the window should be taken 
as the width of the energy distribution SE. Indeed, as long 
as the width squared is extensive (as required in most equi- 
librium situations by the central limit theorem) the difference 
between S cq for any thermodynamic ensemble and S m is not 
only subextensive but actually a constant. This choice thus 
significantly reduces finite size corrections to the entropy. 



Basically, all applications of the entropy to thermodynam- 
ics rely on the fundamental thermodynamic relation: 



dE = TdS — Tdx, 



(9) 



where x is some external parameter like volume or magnetic 
field and T is the corresponding generalized force defined as 
the adiabatic response of the energy to changing such param- 
eter. This relation was first empirically established for quasi- 
static processes in open systems, where TdS describes the 
external heat, and motivated Clausius to introduce entropy in 
the first place. In the context of statistical physics, this relation 
reflects the fact that 



dS 



dE + Tdx 
T 



(10) 



is the proper differential. In other words, the R.H.S. of 
Eq. ( [Tol l does not depend on the details of the process which 
leads to changes of energy and external parameter. For large 
systems, by T one understands the microcanonical tempera- 
ture: T- 1 = dRh\(tt(E)), where tt{E) is the many-body 
density of states QJJ] . For small isolated systems coupled to a 
thermal bath, the relation ( [Tot works if the process connects 
two Gibbs ensembles close in energy and in external parame- 
ter J2l- Then by T in Eq. dTOb one understands the canonical 
temperature. Using the concept of diagonal entropy, it is easy 
to show that in order for Eq. (fTOl to hold for arbitrarily small 
systems it is sufficient that only the initial ensemble be de- 
scribed by the Gibbs distribution JJ]. 

The results of our work reported in the main part of the 
paper @] allow us to establish a very important result that 
in large thermally isolated systems, where the difference be- 
tween canonical and microcanonical statistical ensembles is 
not important, thermalization understood in terms of relax- 
ation of observables to thermal values also implies the validity 
of the fundamental relation (|9). We established that if the inte- 
grability breaking perturbation is sufficiently big then, follow- 
ing the quench, the microscopically defined diagonal entropy 
coincides with the thermodynamic entropy, which is in turn a 
unique function of energy. This means that the time averaged 
state of a closed system after a quench behaves like a thermal 
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state not only from the point of view of statistical physics but 
also from the point of view of thermodynamics. It is interest- 
ing that the onset of the eigenstate thermalization hypothesis 
and the onset of validity of the fundamental thermodynamic 
relation in terms of the strength of the integrability breaking 
perturbation is the same (compare Fig. 1 in the main text and 
Fig. 2 in Ref. H). 

Actually, using the expression of the diagonal entropy 
through the energy distribution we can find leading finite size 
corrections to the fundamental relation. As we established in 
the main paper in ergodic systems only the smooth part of 
the energy distribution contributes the entropy. In the next sec- 
tion we show that even for relatively small systems of size N 
the energy distribution can be well approximated by a Gaus- 
sian and thus the expression in Eq. (|8]l can be applied to the 
relaxed state after the quench. To find the leading 1/7V cor- 
rection to the fundamental relation we imagine that starting 
from this relaxed state the system undergoes a process where 
the average energy E, the external parameter x, and the width 
SE 2 = a 2 undergo some small changes. Note that because 
the system is closed the width of the distribution is in general 
non-canonical (see Refs. 0,01 for additional discussion). 

Expanding the expression for the diagonal entropy in the 
Gaussian approximation [Eq. ([H}] to the first order differen- 
tials in SE, Sa and Sx we find 



SS d ps (3SQ 



1 



1 



Sa 2 



(11) 



where SQ = SE + FSx is the non-adiabatic part of the energy 
change or simply heat ^ and T = 1//3 is the microcanonical 
temperature corresponding to the mean energy in the system. 
We see that if the energy distribution has a canonical width, 
a 2 = a 2 , then the last term in the equation above vanishes 
and the fundamental relation works up to the 1/N 2 correc- 
tions. If the width of the distribution is non-canonical then the 
fundamental relation clearly has a sub-extensive correction, 
which might be significant in mesoscopic size systems. This 
correction is proportional to the change of the width of the 
distribution, which in general is unrelated to SQ. However, if 
we are dealing with thermally isolated ergodic systems then 
Sa 2 and SQ are not independent. They are connected by the 
fluctuation-dissipation relation, which is a direct consequence 
of Jarzynski and Crook's fluctuation relations extended to ar- 
bitrary non-canonical systems (see Ref. H for details). In 
particular, under fairly generic conditions of not very wide 
work distribution and to the order 1/N, we have 

Sa 2 ps 2TSQ. (12) 

Substituting this expression into Eq. (fTTT i we find that the fun- 
damental relation holds for arbitrary distribution with a renor- 
malized temperature: 

T eS SS d = SE + TSx = SQ, (13) 



where 



r n1 



T 



l + T 2 (a- 2 



(14) 



In Ref. ]7[] it was shown that under conditions of validity of 
fluctuation dissipation relation the energy distribution of a re- 
peatedly driven systems acquires universal form with the ratio 
a 2 /a 2 q ps a, where a is a universal constant of the order of 
unity (in some situation this constant can even diverge). Then 
using that a 2 q = T 2 C V we find that the temperature correction 
becomes universal: 



T 



cir 



1-a 1 

a C v 



(15) 



This correction becomes more significant at lower tempera- 
tures where the specific heat is small. 

The relation (fT3l has direct experimental consequences. 
In particular, it implies that the ratio 5Q/T e g = (SE + 
J-5x)/T e B is independent of the details of the dynamical pro- 
tocol, which is used to induce changes in energy and external 
parameter. Likewise using the concept of diagonal entropy 
and the fundamental relation one can make other testable pre- 
dictions in finite size ergodic systems which are not in the 
Gibbs state. In particular, because the effective tempera- 
ture T e g enters the fundamental relation, it should also en- 
ter the Laplace transform needed in order to define free en- 
ergy: F = E — T e gSd- For the same reason, the effective 
temperature will enter all thermodynamic relations including 
the Maxwell's relations. Similarly using the concept of di- 
agonal entropy one can analyze corrections to the thermody- 
namic relations when the system is weakly nonintegrable and 
the stochastic part of the energy distribution gives significant 
corrections to the entropy. This analysis will be the focus of a 
future work. 



ENERGY DISTRIBUTION 

Figure 2 in the main text showed that the distribution func- 
tion of energy for HCBs approaches a Gaussian shape as the 
system transitions to the chaotic limit. The values of the fit- 
ting parameters a, b, as well as E\ n i, and SE used in that figure 
for t', V = 0.00, 0.02, 0.32, 0.48 are given in Table U and HU 



TABLE I: Parameters of Fis. 2 for t u 



0.5, V m 



2.0 (top panels) 



t' = V 


a 


b 


Eim 


SE 


0.00 


-4.43 


0.92 


-3.84 


1.60 


0.02 


-4.40 


0.92 


-3.86 


1.54 


0.32 


-4.80 


1.76 


-4.62 


1.87 


0.48 


-5.57 


1.77 


-5.33 


1.81 



' cq 



In the figure below, we consider again the case of HCBs 
and show both the diagonal elements of the density matrix, 
p nn , vs energy (top panels) and semi-logarithmic plots of the 
distribution functions of energy (bottom panels). The quench 
is from i; n j = 0.5, Vim = 2.0 and two values of t',V' are 
considered: in the integrable regime, t' = V' = 0.00 (left 
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TABLE II: Parameters of Fig. 2 for tu = 2.0, Vm = 0.5 (bottom 
panels) 



f = V 


a 


6 


Em 


SE 


0.00 


-4.00 


0.06 


-3.84 


0.70 


0.02 


-4.00 


0.06 


-3.86 


0.71 


0.32 


-4.66 


0.62 


-4.62 


0.78 


0.48 


-5.34 


0.75 


-5.33 


0.88 



panels), and in the chaotic domain, t' = V = 0.32 (right 
panels). 

As seen from the top panels, there are large fluctuations 
of p nn and no clear features distinguishing the integrable and 
chaotic regimes, except probably that larger fluctuations are 
noticeable in the integrable case. Correlations between those 
fluctuations and the largely fluctuating expectation values of 
observables in the eigenstates of the Hamiltonian preclude 
thermalization in integrable systems (see Ref. [9] in the main 
text). 

The semi-logarithmic plots in the bottom panels emphasize 
that in the nonintegrable case, not only the center of W(E) 
approaches a Gaussian shape, but the behavior of the tails of 
the distribution is also Gaussian (see Ref. [14,16] in the main 
text). 




E 



FIG. 4: (Color online) Weights of the components of the initial state 
(top panels) and normalized distribution function of energy (bottom 
panels). Bosonic system, L — 24, T = 3.0, and quench from t-mi — 
0.5, Mm = 2.0. Left panels: integrable domain (t' = V' = 0.00); 
right panels: chaotic regime (t' — V' — 0.32). Solid smooth line: 
best Gaussian fit (y^^a) -1 ^^ -6 ' // ' 2a ' for parameters a and 6; 
dashed line: (y/^SE^e-^-^ 2 /(2Se2) . The values of a, b, 
E m i and 5E are given in Table U 
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